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Main  contributions  of  the  Research  Project 

Consider  a  rocket,  modeled  by  a  cylinder,  rotating  and  moving  at  high  speed  strikes 
a  target.  It  generate  a  longitudinal  stress  a  and  torsional  stress  r. 


If  the  material  is  a  generally  nonlinear  isotropic  elastic  solid,  the  governing  equations  are  a 
system  of  hyperbolic  conservation  laws 

Uj  -  P(U),  =  0  (1) 

which  also  arises  in  other  physical  problems  such  as  in  fluid  dynamics  and  oil  recovery 
problems.  The  function  form  for  F  may  be  different  for  other  physical  problems  but  the 
underlined  mathematical  problems  are  the  same.  Thus  the  results  obtained  from  this 
project  and  from  other  researchers  on  fluid  dynamics  and  oil  recovery  complement  each 
other.  The  main  contributions  from  the  present  project  are  as  follows. 

(1)  Classification  of  2x2  non-strictly  hyperbolic  system.  For  fluid  dynamics,  Eq. 
(1)  is  a  2x2  system.  For  waves  in  solids,  it  can  be  reduced  to  a  2x2  system  if  the  boundary 
conditions  are  prescribed  in  terms  of  stress  rather  on  velocity.  There  are  two  wave  speeds 
Cj  and  Cg  and,  as  it  happens  in  many  applications,  c^^  =  Cg  for  certain  U  which  is  called 
the  "umbilic"  point.  The  system  is  then  "non-strictly  hyperbolic".  Non— strictly 
hyperbolic  systems  have  been  the  subject  of  intensive  study  by  many  applied 
mathematicians  and  fluid  dynamidsts.  One  of  the  fundamental  problems  for  non— strictly 


□  □ 
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hyperbolic  systems  is  the  geometry  of  the  "wave  curves"  near  an  umbilic  point.  The  2x2 
system  can  be  classified  iccorriiv.g  to  the  geometry  of  the  wave  curves  near  the  umbilic 
point.  For  example,  the  wave  curves  for  nonlinear  isotropic  elastic  solids  whose  strain 
energy  is  a  function  of  (a,  r)  of  order  up  to  three,  may  have  the  geometry  shown  in  Fig.  1 


Point  0  is  the  zero  stress  state.  If  the  stress  is  initially  at  point  B’  and  the  boundary 
condition  at  x  =  0  is  at  point  A,  the  stress  variation  at  a  fix  location  x  is,  as  time  t 
increases,  follows  the  wave  curve  B’MA.  The  corresponding  wave  pattern  in  the  (x,  t) 


Fig.  la 


Fig.  lb 


plane  consists  of  two  simple  wave  fans  as  shown  in  Fig  la.  In  Fig.  1,  the  solid  lines  denote 

wave  curves  associated  with  the  fast  wave  speed  c^^  and  tht  aotted  lines  with  the  slow  wave 

speed  Cg.  One  can  only  move  along  the  direction  of  the  arrows  and  can  change  from  a  solid 

line  to  a  dotted  lines.  If  the  path  calls  for  changing  from  a  dotted  curve  to  a  solid  curve  or 

going  against  the  direction  of  the  arrows,  a  shock  wave  will  be  generated.  For  instance,  if 

the  initial  condition  is  at  point  B  and  the  boundary  condition  is  at  A  a  shock  wave  is 

generates  as  shown  in  Fig.  lb.  The  region  before  the  shock  wave  (denoted  by  a  double  line) 

is  the  constant  stress  state  B  and  the  solution  behind  the  shock  wave  is  at  A. 

* 

The  point  a  in  Fig  1  is  the  umbilic  point  at  which  c^  =  c^.  Depending  on  the  the 
material  constants,  the  wave  curves  have  different  geometry  as  shown  in  Figs.  2—4.  Figs. 
1—5  are  based  on  the  2x2  hyperbolic  system  of  conservation  laws.  The  classification  of  the 
geometry  of  the  wave  curves  is  one  of  the  main  contributions  of  this  research.  If  the  strain 
energy  is  assumed  to  depend  on  (<r,r)  of  order  up  to  four,  the  umbilic  point  may  not  be  at 
the  a  axis.  The  wave  curves  near  the  umbilic  point  depend  on  two  composite  material 
parameters  m  and  k.  Figures  1—4  are  special  cases  in  which  m  =  0.  The  geometry  of  the 
wave  curves  near  the  umbilic  point  can  be  classified  into  five  cases  as  shown  in  Fig.  5. 


Fig.  2.  Simple  wave  curves  for  Case  2:  - 1  <  it  s  0 


Fig.  3  Simple  wave  curves  for  Case  3:  0<i  s  I  For  Case  3a.  m,  >mj>o*  For  Case  3b, 

/n,  <  /n,  <  o'. 


'Hie  Reiinann  probletn  of  pUae  waves 


Fig.  S  Classification  according  to  the  number  of  wave  lines  (WL)  and  inflection  lines  (IL). 


This  classification  is  simpler  than  what  is  available  in  the  literature.  We  did  more.  We 
classify  in  more  detail  witLin  each  case.  For  instance,  Case  2  in  Fig.  5  can  be  further 
classified  as  shown  in  Fig  6. 
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Fig.  Further  classification  of  Case  2. 


The  geometry  of  wave  curves  in  each  subdivision  is  of  course  different.  The  implication  is 
that,  in  solving  the  Riemann  problem,  it  would  eliminate  the  problem  of  guessing  which 
wave  curve  geometry  to  choose  for  the  solution. 

(2).  Wave  curves  with  the  presence  of  an  umbilic  point  and  an  umbilic  line.  For  the 

2x2  system  with  the  assumption  that  the  potential  function  is  of  the  third  order  in  U, 

there  is  only  one  umbilic  point  but  no  umbilic  line.  In  applications  one  may  want  to 

include  fourth  and  higher  order  terms  for  a  better  accuracy.  The  system  then  may  have 

more  than  one  umbilic  point  and/or  an  umbilic  line.  For  instance,  with  the  fourth  order 

term  Fig.  1  is  replaced  by  Fig.  7.  Point  O  in  Fig.  7  is  the  umbilic  point  and  hence  is 

* 

identical  to  the  umbilic  point  a  in  Fig.  1.  The  wave  curves  near  point  0  in  Fig.  7  are 

* 

identical  to  that  near  point  a  in  Fig.  1.  However  wave  curves  away  from  the  umbilic 
point  in  the  two  figures  are  different.  We  have  studied  several  potential  functions  of  fourth 
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order  in  U  and  indeed  found  many  interesting  phenomena.  For  instance,  in  Fig.  7,  there  is 
an  umbilic  line  and  also  an  inflection  line  across  which  the  wave  speed  attains  a  maximum 
or  a  minimum.  The  corresponding  problems  in  fluid  dynamics  do  not  seem  to  have  been 
studied.  There  is  no  doubt  that  some  of  the  results  obtained  are  applicable  to  other 
physical  problems. 

(3)  Generalized  Riemann  problem.  The  standard  Riemann  problem  for  waves  in 
solids  is  that  the  initial  conditions  at  t  =  0,  x  >  0  are  constant  and  the  boundary 
conditions  at  x  =  0,  t  >  0  are  another  constant.  This  is  sometimes  call  the 
Goursat— Riemann  problem.  If  the  stress  (a,  t)  at  x  =  0  after  the  impact  is  not  a  constant 
(which  is  the  case  in  most  practical  problems)  we  may  assume  that,  for  small  t, 

<t(0,  t)  =  (t(0,  0)  +  a  t, 
r(0,  t)  =  r(0,  0)  +  ft, 

where  f,  f  are  constant.  This  is  a  generalized  Riemann  problem.  It  is  a  first  order 
approximation  in  the  boundary  conditions.  We  could  also  consider  second  and  higher  order 
approximations  but  for  the  near  field  solutions  near  the  impact  point  x  =  t  =  0,  the  first 
order  approximation  is  suflident.  The  generalized  Riemann  problems  have  only  recently 
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attracted  the  attentions  of  French  applied  mathematicians.  As  it  turns  out,  our 
approaches  in  solving  the  generalized  Riemann  problem  which  are  based  ci  my  earlier 
investigations  on  plastic  waves  about  ten  years  ago,  coincide  with  the  approaches  suggested 
by  the  French  applied  mathematicians.  However,  they  did  not  apply  their  approaches  to 
any  physical  problems.  We  have  been  studying  generalized  Riemann  problem  with 
applications  to  impact  on  solids  in  the  past  ten  months  or  so.  Preliminary  results  have 
been  reported  at  a  workshop  on  shock  waves  and  are  appearing  in  one  of  the  volumes  in 
SIAM  Proceedings  on  Applied  Mathematics.  More  results  are  being  compiled  and  will  be 
reported  in  the  near  future.  To  show  the  difference  in  the  solution  for  the  Riemann 
problem  and  for  the  generalized  Riemann  problem,  consider  the  wave  curve  BA  in  Fig.  1 
which  results  in  a  shock  wave  as  shown  in  Fig.  la.  For  the  generalized  Riemann  problem, 
the  solution  is  shown  in  Figs.  8a  or  8b  depending  on  the  values  of  a,  f.  In  Fig.  8a  there  is  a 
non-steady  shock  wave  with  increasing  speed  and  amplitude  followed  by  a  weak 
discontinuity  along  c^  characteristics.  In  Fig.  8b  the  shock  wave  is  not  affected  but  the 
region  after  c^  is  not  a  constant  region.  It  is  a  simple  wave  region. 


Fi;^.  3a 


Fig.  8b 


9 


Publications  under  this  project. 


[1]  Guangshan  Zhu  and  T.  C.  T.  Ting,  "Classification  of  2x2  non-strictly  hyperbolic 
systems  for  plane  waves  in  isotropic  elastic  solids,"  Int.  J.  Eng.  Sci.  27,  1621—1638 
(1989). 

[2]  T.  C.  T.  Ting  and  Tankin  Wang,  "The  Goursat— Riemann  problem  for  plane  waves 
in  isotropic  elastic  solids  with  velocity  boundary  conditions,"  in  Multidimensional 
Hyperbolic  Problems  and  Computations.  The  IMA  Volumes  in  Mathematics  and 
Its  Applications,  29.  James  Glimm  and  Andrew  J.  Majda,  eds.  Springer— Verlag. 
367-386  (1991). 

[3]  Guangshan  Zhu  and  T.  C.  T.  Ting,  "Wave  curves  with  the  presence  of  an  umbilic 
point  and  an  umbilic  line  for  plane  waves  in  isotropic  elastic  solids,"  Wave  Motion. 
13,  277-290  (1991). 

[4]  Guangshan  Zhu  and  T.  C.  T.  Ting,  "Hyperbolic  conservation  laws  with 
self— similar  wave  curves,"  Proc.  Third  Int.  Conference  on  Hyperbolic  Problems, 
Theory,  Numerical  Methods  and  Applications,  (held  in  Uppsala,  Sweden,  June 
11-15,  1990.)  B.  Engquist  and  B.  Gustafsson,  eds.  Studentlitteratur,  Lund, 
Sweden,  1010-1024  (1991). 

[5]  Tankin  Wang  and  T.  C.  T.  Ting,  "Growth  or  decay  of  shock  waves  in  the 
generalized  Goursat-Riemann  problem,"  in  Viscous  Profiles  and  Numerical 
Methods  for  Shock  Waves.  M.  Shearer,  ed.  SIAM  Proceedings  Series.  SIAM, 
Philadelphia.  161-174  (1991). 

[6]  Tankin  Wang  and  T.  C.  T.  Ting,  "The  generalized  Goursat-Riemann  problem  for 
plane  waves  in  isotropic  elastic  solids,"  under  preparation. 


Copies  of  publications  [1]  and  [5]  are  attached  herewith. 


Int.  J.  Engitt  Sa.  Vol.  27.  No.  12.  pp.  1621-1638.  1989 
Primed  in  GftM  Briuin.  All  righa  reserved 


0020-7225/89  S3.00  +  0.00 
Copyright  ©  1989  Pergamon  Press  pic 


CLASSmCATION  OF  2  x  2  NON-STRICTLY  HYPERBOLIC 
SYSTEMS  FOR  PLANE  WAVES  IN  ISOTROPIC  ELASTIC 

SOLIDS 

GUANGSHAN  ZHU  and  T.  C.  T.  TING 

Depanmeot  of  Civil  Engineering,  Mechanics  and  Metallurgy,  University  of  Illinois  at  Chicago.  Box 

4348.  Chicago.  IL  60680.  U.S.A 

Abetract — The  governing  equations  for  plane  waves  in  generally  nonlinear  isotropic  elastic  solids  are 
a  system  of  6  x  6  hyperbolic  conversation  laws.  For  the  half-space  Riemann  problem  in  which  the 
initial  conditions  at  r  =  0,  x  >  0  and  the  boundary  conditions  at  x  =  0.  r  >  0  are  constant,  the  system  is 
equivalent  to  3  x  3  system  in  the  full-space  Riemann  problem.  It  is  further  reduced  to  a  2  x  2  system 
due  to  the  fact  that  one  of  the  characteristic  wave  speeds  is  linearly  degenerate.  For  hyperelastic 
materials  for  which  there  exists  a  potential  W  whose  gradients  provide  the  strains,  the  wave  curves 
near  an  isolated  umbilic  point  are  represented  by  the  potential  of  the  form 

Wfa,  t)  =  +  —^1^ 

which  contains  two  parameters  k  and  m.  The  classification  of  the  geometry  of  wave  curves  depends  on 
the  values  of  k  and  m  and  can  be  classified  into  five  cases.  The  potential  function  considered  here  is 
equivalent  to 


IV  =  I  +  ba^x  +  or 

considered  by  Schaeffer  and  Shearer  where  a  and  b  are  parameters.  The  classification  presented  here 
seems  to  provide  simpler  algebraic  expressions.  It  also  renders  more  refined  classification  as  shown  in 
the  paper. 


1.  INTRODUCTION 


Consider  a  plane  wave  propagating  in  the  Xi-direction  in  a  fixed  rectangular  coordinate  system 
X,,  xj,  X3.  The  material  occupies  the  half-space  XiarO.  Let  o,  Tj,  T3  be,  respectively,  the 
normal  stress  and  the  two  shear  stresses  on  a  plane  Xi  =  constant.  If  u,  Vj,  V3  are,  respectively, 
the  particle  velocity  in  the  x^,  X2,  X3  direction,  the  equations  of  motion  and  the  continuity  of 
displacements  can  be  written  as  a  6  x  6  system  of  hyperbolic  conservation  laws  [1,2] 


U.-F{U),  =  0. 

U  =  (a,  T2,  Tj,  u.  V2,  Vj). 

F(U)  =  (pu,  pv2,  pv„  e,  72,  73).^ 


(1.1) 


In  the  above,  x  =x,,  t  is  the  time,  p  is  the  mass  density  in  the  undeformed  state,  and  s,  72,  7, 
are,  respectively,  the  extensional  strain  in  the  Xi-direction  and  the  two  shear  strains.  The 
constitutive  equations  of  the  material  specify  the  relations  between  a,  T2,  Tj  and  e,  Yi,  Yj-  For 
general  nonlinear  elastic  materials  for  which  a,  Tz,  t,  are  known  functions  of  e,  73,  73,  the 
system  of  equations  (1.1)  has  six  characteristic  wave  speeds  ±Cj,  i  =  1,  2,  3  [3,4].  Without  loss 
in  generality,  we  let 


0<  C3  <C2  2  c,. 


The  Riemann  problem  is  a  special  Cauchy  problem  in  which  the  initial  conditions  at  ( =  0. 
x>0  and  at  f  =  0.  x<0  are  constant  [1,5-9].  For  plane  waves  in  a  half-space,  the  Riemann 
problem  prescribes  constant  initial  conditions  at  r  =  0.  x  >  0  and  constant  boundary  conditions 
at  X  =0,  f  >0.  The  solution  in  the  (x,  t)  plane  depends  on  one  parameter  x/f  only.  If  the 
solution  is  continuous  in  xlt,  we  have  a  simple  wave  solution  in  which  x/f  =  c,  the  characteristic 
wave  speed.  In  particular,  the  stresses  a,  T2,  T3  are  continuous  functions  of  c  and,  as  c  varies, 
the  stresses  trace  out  a  “simple  wave  curve”  in  the  stress  space.  Since  we  are  considering  the 
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half-space  Ji^O,  we  need  only  the  positive  wave  speeds  c,,  Cj,  c,.  We  therefore  have  three 
families  of  simple  wave  curves  and  (1.1)  for  the  Riemann  problem  in  the  half-space  is 
equivalent  to  a  3  x  3  hyperbolic  system. 

The  solution  to  the  Riemann  problem  is  of  fundamental  importance  in  finding  the  solution  to 
a  more  general  initial  and  boundarv  vpIup  problem  [10].  The  solution  requires  precise 
classification  of  the  simple  wave  curves.  The  lUemann  problem  appears  not  only  in  wave 
propagation  in  solids,  but  also  in  other  physical  applications  such  as  oil  recovery  problems 
[10.11]  as  well  as  in  solving  two-dimensional  Riemann  problems  [12.13].  In  finding  the  solution 
to  the  Riemann  problem  we  need  the  knowledge  of  the  geometry  of  simple  wave  curves.  This 
leads  to  the  classification  of  the  hyperbolic  system  according  to  the  geometry  of  simple  wave 
curves.  In  the  rest  of  the  paper,  we  refer  to  simple  wave  curves  as  wave  curves  for  brevity 
because  we  will  not  be  concerned  with  shock  wave  curves. 

In  Section  2  we  point  out  that  when  the  material  is  isotropic,  (1.1)  is  linearly  degenerate  with 
respect  to  C2-  It  suffices  therefore  to  consider  wave  curves  associated  with  c,  and  only  and 
(1.1)  is  reduced  to  a  2x2  system.  For  hyperelastic  materials  for  which  the  strains  are 
obtainable  from  the  gradients  of  potential  function  W,  the  wave  curves  associated  with  c,  and 
C3  are  orthogonal  to  each  other.  This  is  presented  in  Section  3.  The  classification  of  the  2  x  2 
system  is  studied  in  Sections  4-7.  For  the  system  under  consideration,  there  exists  an  umbilic 
point  at  which  Cj  =  Cj.  The  system  therefore  is  not  strictly  hyperbolic  [14-18].  For  second  order 
hyperelastic  materials  for  which  W  can  be  approximated  by  a  polynomial  in  stress  of  order  up 
to  three,  there  is  one  umbilic  point  which  is  located  on  the  a-axis.  The  classification  of  2  x  2 
systems  for  second  order  materials  is  presented  in  Section  4.  In  Section  5  we  consider  the  third 
order  materials.  There  are  now  more  than  one  umbilic  point.  We  may  also  have  an  umbilic 
line.  If  the  umbilic  point  is  on  the  c7-axis,  the  classification  is  identical  to  that  for  the  second 
order  materials.  If  the  umbilic  point  is  not  on  the  a-axis,  we  show  in  Section  6  that  the 
classification  is  identical  to  that  of  2  x  2  systems  of  general  potential  function  near  an  isolated 
umbilic  point.  Section  7  contains  the  main  results  of  our  paper.  We  present  the  classification  of 
2x2  systems  of  general  potential  function  near  an  isolated  umbilic  point.  The  classification  is 
different  from  that  in  [11]  but  can  be  shown  to  be  equivalent  if  one  rotates  the  coordinate  axes 
of  the  stress  space.  The  classification  presented  here  is  more  refined  and  appears  to  offer 
simpler  algebraic  expressions  for  the  boundaries  between  different  cases. 


2.  EQUIVALENT  TO  2  x  2  SYSTEMS  FOR  ISOTROPIC  ELASTIC  MATERIALS 


When  the  material  is  isotropic,  the  relation 


T2/r3  =  Yi/y., 


holds  [8].  If  we  let 

T;  =  T  cos  6,  T3  =  T  sin  B,  t  >  0. 

72  =  y  cos  6,  73  =  y  sin  6,  y  >  0, 

the  rectangular  coordinate  system  (a,  r,,  r,)  is  transformed  into  a  cylindrical  coordinate  system 
(a,  r,  6).  It  can  be  shown  [8,9]  that  the  wave  curves  associated  with  c,  and  C3  are  plane 
polarized  on  the  6  =  constant  plane.  On  the  other  hand,  the  wave  curves  associated  with  c,  is 
circularly  polarized  on  the  circle 


a  =  constant,  t  =  constant,  for  c  =  C2. 


Moreover,  C2  is  a  constant  on  the  circular  wave  curve.  Therefore  C2  is  linearly  degenerate  [10] 
and  the  simple  wave  is  in  fact  a  shock  wave.  For  the  Riemann  problem  in  which  the  initial 
condition  and  the  boundary  condition  for  6  are  different,  all  we  have  to  do  is  to  introduce  a  C2 
shock  wave  which  changes  the  value  of  6  from  the  initial  value  to  the  value  prescribed  at  the 
boundary.  This  is  one  of  few  special  cases  in  nonlinear  problems  in  which  the  principle  of 
superposition  applies. 
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It  is  therefore  sufficient  to  consider  6  to  be  constant  and  the  6  x  6  system  of  (1.1)  is  reduced 
to  a  4  X  4  system 

U,-F(lJ),  =  0. 


U  =  (a,  T,  u,  u), 

F{U)~(pu,  pit.  E,  y).  J 
The  characteristic  wave  speeds  c,,  c,  are  the  eigenvalues  of  the  eigenrelation 

(I  +  cA)U'  =  0. 


(2.1) 


oJ-  «-u  jj 


(2.2) 


In  the  above,  I  is  a  unit  matrix,  a  prime  denotes  differentiation  with  c  and  the  subscripts  o  and 
T  denote  differentiation  with  these  variables.  Here  we  have  assumed  that  the  relations  between 
a,  T  and  e,  y  are  invertable.  If  not,  we  limit  our  attention  to  the  region  where  they  are 
invertable. 

Since  G  depends  on  a  and  r  only,  the  four  equations  in  (2.2)i  can  be  reduced  to  two 
equations  as 


(G  -  ql)s'  =  0.  ] 

s=(a.  T),  ij  =  (pc’)-'.  ' 
If  we  consider  the  2  x  2  hyperbolic  system 

U,  +  F(U).  =  0. 

U  =  (a.  r).  F(U)  =  (£.  y), 


(2.3) 


(2.4) 


the  Riemann  problem  for  this  system  also  leads  to  (2.3).  Therefore,  the  Riemann  problem  for 
the  half-space  of  4  x  4  system  (2.1)  is  equivalent  to  the  Riemann  problem  for  the  full-space  of 
2x2  system  (2.4). 


3.  HYPERELASTIC  MATERIALS 

For  hyperelastic  materials  there  exists  a  potential  W,  the  complementary  strain  energy  (19), 


such  that 

y  =  W, 

The  matrix  G  in  (2.2)  has  the  form 

G  =  [ 

LW„, 

w„y 

which  is  symmetric.  The  eigenvalues  r;  of  G  are 

h.=^{(W„„  +  W'„)-y}=(pcf)-‘, 

Vy  =  \{(Ka  +  W„)+Y}=(pcl)-\ 

The  eigenvector  s’  of  (2.3)  provides  the  tangent  of  wave  curves  in  the  stress  space  (o  r).  The 
differential  equation  for  wave  curves  is  therefore 

dr  2W„,  -(w„„-W„)Ty 

0’  do  {W„„-W„)TY  2W„, 


(3.3) 


1624 


GUANGSHAN  ZHU  and  T.  C.  T.  TING 


in  which  the  upper  sign  (or  the  lower  sign)  is  for  the  c,  (or  c^)  wave  curves.  Due  to  the 
symmetry  of  G,  the  eigenvectors  associated  with  Ci  and  are  orthogonal  to  each  other.  This 
means  that  the  wave  curves  for  Cj  and  Cj  are  orthogonal  to  each  other.  Regarding  z  =  W(a,  r) 
as  a  surface,  it  can  be  shown  [11]  that  the  wave  curves  are  the  projection  of  the  lines  of 
curvature  of  the  surface  on  the  (a,  t)  plane. 

When  Cl  =  C3,  i.e.,  rji  =  r/j,  the  system  is  not  strictly  hyperbolic.  The  sfess  (a,  v)  at  which 
Cl  =  C3  is  called  the  umbilic  point.  We  see  from  (3.2)  that  the  umbilic  point  is  located  at 

=  w„  =  0.  (3.4) 

If  (3.4)  are  satisfied  for  a  one-parameter  family  of  points,  we  have  an  umbilic  line.  Since  the 
umbilic  line  is  a  line  of  curvature  of  the  surface  z  =  W(a,  r),  the  umbilic  line  is  a  wave  curve. 
We  also  see  that  at  the  umbilic  point  or  the  umbilic  line,  the  wave  curves  for  Ci  and  C3  may  not 
be  orthogonal  to  each  other  [8,9]. 

The  differential  equation  (3.3)  for  the  wave  curves  has  the  following  invariance  properties: 

(i)  If  we  multiply  or  divide  W  by  a  constant,  (3.3)  remains  unchanged. 

(ii)  Equation  (3.3)  is  invariant  with  a  translation  and/or  a  rotation  of  the  coordinate  axes 
a  and  t  [11]. 

(iii)  If  we  expand  W{a,  t)  in  powers  of  a  and  r,  (3.3)  remains  the  same  if  we  ignore  the 

constant  terms  and  the  linear  terms.  It  also  remains  the  same  if  we  ignore  the  term 
(tr-t-  T^)  because  this  term  contributes  nothing  to  and  which  appear  in 

(3.3). 

The  above  invariance  properties  form  the  bases  in  the  sequel  for  reducing  the  potential 
function  W  to  an  equivalent  form  for  classification. 


4.  CLASSIFICATION  FOR  SECOND  ORDER  HYPERELASTIC  MATERIALS 


From  (3.1),  if  the  strains  are  to  be  expressed  in  powers  of  a,  t  of  order  up  to  two,  we  have  to 
expand  W  in  terms  of  a  and  x  of  order  up  to  three.  We  may  ignore  the  constant  term  in  W 
since  it  contributes  nothing  to  the  strains  c  and  y.  We  ignore  also  the  first  order  terms  in  a  and 
T  because  the  strains  would  otherwise  be  non-zero  when  the  stresses  are  zero.  Noticing  that  for 
isotropic  materials  W  should  be  an  even  function  in  the  shear  stress  r,  we  have,  for  the  second 
order  hyperelastic  materials, 

a  d  b  ,  € 

W(a,  t)  =  - 0-  +  - T^  +  -cf^  + -or,  (4.1) 

in  which  a,  d,  b,  e  are  constants.  From  (3.1), 


e  =  ao  +  -  (ba^  +  ex^), 


(4.2) 


Y  =  dx  +  eax,  J 

and  hence  a  and  d  are  the  elastic  constants  for  linear  materials  while  b  and  e  are  the  second 
order  elastic  constants.  The  constants  a  and  d  are  related  to  Lame  constants  A  and  fi  by 

a~'  =  k  +  2fi,  d~' -  n. 

With  W  given  by  (4.1),  the  umbilic  point  obtained  from  (3.4)  is  located  at 

0=0.,  T  =  0. 


where 


o. 


d  -  a 
b  -  e 


Therefore,  the  umbilic  point  is  on  the  0-axis. 
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Equation  (4.2)  tells  us  that,  when  e  =  0,  the  stress-strain  laws  for  tension  and  shear  are 
uncoupled.  The  2x2  system  is  then  decoupled  into  two  1x1  systems.  This  special  case  has 
been  studied  in  [20].  We  assume  here  that  e  #0.  Next,  we  see  that  if  we  change  the  sign  of  e. 
(4.1)  and  (4.2)  remain  the  same  if  we  also  change  the  sign  oi  b,  a  and  c.  The  wave  curves  for 
e  <  0  can  therefore  be  obtained  from  that  for  e  >  0  by  reversing  the  direction  of  the  a-axis  and 
the  sign  of  b.  Hence  it  suffices  to  consider  e  >0.  Finally,  using  the  invariance  properties  (i),  (ii) 
and  (iii)  discussed  in  Section  4  we  divide  W  by  c  and  move  the  origin  of  (a,  t)  plane  to  (a,,  0). 
Equation  (4.1)  is  then  replaced  by 

1  _  ^  1 

b 

k=l  — ,  a  =  o~  a», 

e 

and  the  differential  equation  (3.3)  for  the  wave  curves  can  be  written  as 

-<PT(1  +  <PY'"={4>T(1  +  4>")*'^}->.  (4.4) 


which  contains  only  one  parameter  k. 

Equation  (4.4)  can  be  integrated  in  a  closed  form  [21].  The  wave  curves  have  different 
geometry  depending  on  the  value  of  k.  They  can  be  classified  into  four  cases  as  follows  [8,  9]: 

Case  1:  k<  -1. 

Case  2:  —  1  s  ifc  <  0. 

(4.5) 

Case  3;  0  <  A:  s  1. 

Case  4:  1  <  Ac. 

The  geometry  of  wave  curves  for  each  case  and  the  solution  to  the  Riemann  problem  have 
been  discussed  in  detail  in  [8,9].  We  will  see  that  the  classification  (4.5)  is  a  special  case  of  the 
classification  of  general  2x2  systems  to  be  considered  in  Section  7. 


5.  CLASSIFICATION  FOR  THIRD  ORDER  HYPERELASTIC  MATERIALS 

We  now  include  the  third  order  terms  in  a  and  x  for  the  strains  e  and  y.  This  means  that  W 
must  contain  the  fourth  order  terms  and  we  write 

W  =  ^(7^  +  ^T^  +  ^(T^  +  ^ar^  +  ^6,a-‘  +  ^  ^  (5. 1 ) 

in  which  du  bi,  63  are  the  third  order  elastic  constants.  The  strain-stress  relations  (3.1)  have 
the  expressions 

e  =  ao  +  ^  (bo^  ex^)  +  6,<r^  -I-  d^x'^a, 

y  =  dx  +  eax  +  Q  63^  +  630^^  r. 

There  are  now  more  than  one  umbilic  point  and  there  may  be  umbilic  lines.  With  (5.1).  the 
conditions  (3.4)  for  the  umbilic  point  are 

a  +  bo  +  b\C^  +  bjx^  =  d  +  ea  +  bjcr  -P 

r(e  +  2630)  =  0. 
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There  are  several  solutions  to  the  two  equations.  We  discuss  the  solutions  separately  below. 

(A)  Let  at*',  a*,^'  satisfy  the  equation 

a  +  ba  +  6f(r  =  d  +  eo  +  6^cr.  (5.3) 

’Vc  then  have  two  umbilic  points  at  (at",  0)  and  (at^*,  0).  They  are  on  the  a-axis.  To  study  the 
geometry  of  wave  curves  near  the  umbilic  points,  we  rewrite  the  right  hand  side  of  (5.1)  in 
terms  of  a  and  r  where 

a  =  a  —  at",  i  =  1.  2. 

Using  the  invariance  properties  (ii)  and  (iii)  and  ignoring  the  fourth  order  terms,  we  replace 
(5.1)  by 

W  =  2  +  -  dr^, 

6  2 

6  =  b  +  2d, at",  e  =  e  +  26, at". 


Assuming  that  e  #0  and  using  the  invariance  property  (i).  we  have 

W=——d^  +  -aT^,  k=\--.  (5.4) 

6  2  e 

This  is  identical  to  (4.3).  The  geometry  of  wave  curves  near  the  umbilic  point  (af’,  0)  or 
(ot^*,  0)  is  therefore  identical  to  that  discussed  in  the  last  section. 

(B)  If  63^0  and  62  #63,  (5.2)  have  the  solution 


Oo=  - 


To  = 


{(d  -a)  +  (e-  b)ao  +  (63  -  6|)a^} 

(63 -<52) 


(5.5) 


If  To  =  0,  the  umbilic  point  is  on  the  a-axis  and  the  discussion  in  Case  (A)  applies.  If  Tq  #  0,  we 
have  an  umbilic  point  at  (oq,  to)  which  is  not  on  the  a-axis.  Let 

O  —  O  (7q,  T  —  X  Tq- 

Using  the  invariance  properties  (ii)  and  (iii)  and  ignoring  the  fourth  order  terms  we  obtain 

W  =  fd'  + +  -r^, 

6  2  6 

g  =  d  +  2d,ao,  6  =  262X0,  e  =  26310. 

Since  e  #  0,  we  may  divide  W  by  e  to  obtain 

W  =  +  +  (5.6) 

6  2  6 

g  b 

m  =  ^,  k  =  l--. 
e  e 

We  will  show  in  the  next  section  that  for  a  general  potential  function  IV  near  an  isolated 
umbilic  point,  W  also  has  the  expression  (5.6).  The  classification  of  2  x  2  systems  with  W  given 
by  (5.6)  is  studied  in  Section  7. 

(C)  If  62  =  63=^0,  and  (5.5)  satisfies  (5.3).  we  have  an  umbilic  line  at 

a„=  -r^.  r  arbitrarv. 

263 

This  case  is  being  studied  separately  in  [22]. 
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(D)  If  ^3  =  e  =  0,  we  have  the  umbilic  line  given  by 

a  +  ba  +  6, <7^  =  d  4- 

If  62 #0.  the  umbilic  line  is  a  parabola,  hyperbola  or  ellipse  depending  on  whether  6,63  is 
zero,  positive  or  negative.  This  case  has  been  studied  in  (20],  The  degenerate  case  63  =  0  leads 
to  two  straight  umbilic  lines. 


6.  THE  GENERAL  POTENTIAL  FUNCTION  NEAR  AN  ISOLATED  UMBILIC 

POINT 

We  now  consider  a  general  function  W  and  assume  that  (Oo,  ro)  is  an  isolated  umbilic  point. 
This  means,  from  (3.4), 

To)  =  W„(<7o,  To), 

WariOo,  ro)  =  0. 

If  we  expand  W  in  Taylor  series  at  the  umbilic  point,  use  the  invariance  properties  (ii)  and  (iii) 
and  ignore  the  fourth  order  terms,  we  obtain 

I?  ,  c  ,  ft  ,  6  , 

W  =  |a^  +  -a-f +  -dT‘  +  -r.  (6.1) 

^  OQf  X  —  T  Tq, 

In  the  above,  g,  e,  h  and  b  are  the  third  order  derivatives  of  W  at  the  umbilic  point  (oq,  Tq). 
Without  loss  in  generality,  we  may  assume  that  A  =  0.  If  A  #  0,  we  apply  the  invariance 
properity  (ii)  by  rotating  the  coordinate  axes  a,  x  such  that  the  term  disappears.  Hence  we 
consider 

W  =  fd^  +  Vf  +  |f^ 

6  2  6 

The  special  case  e  =  0  decouples  the  2  x  2  system  into  two  1  x  1  systems  [20].  We  therefore 
assume  e  ^0.  Using  the  invariance  property  (i)  to  divide  W  by  e,  we  have 

W  =  -^vVt  +  i^r^  (6.2) 

6  2  6  ’ 

g  b 

m=^,  k  =  \ - . 

e  e 


This  is  identical  to  (5.6). 

Equation  (4.3)  is  a  special  case  of  (6.2)  when  m  =  0  and  the  role  of  a  and  f  are  interchanged. 


7.  CLASSIFICATION  OF  GENERAL  2x2  SYSTEMS  NEAR  AN  UMBILIC 

POINT 


We  now  study  the  classification  of  general  2x2  systems  in  which  the  potential  function  W 
near  an  umbilic  point  is  given  by  (6.2).  For  simplicity,  we  omit  the  overbars  and  write 


I.,/  ^  ^  3 

W(a.  T)  =  -cr  +  -(rT  +  — —  r. 

J  z  o 


(7.1) 


The  umbilic  point  is  now  at  the  origin.  If  we  change  the  sign  of  m  and  a.  W  remains  the  same. 
It  suffices  therefore  to  consider 

m  >0, 
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because  the  wave  curves  for  m<0  can  be  obtained  from  that  for  m  > 0  by  reversing  the 
direction  of  the  a-axis.  The  wave  speeds  c,.  cj  given  by  (3.2)  are 

t;,  =  ^  {2ma  +  (2  -  A:)t  -  V}  =  (pc?)”', 

7/3  =  ^{2ma  + (2  - /c)r  +  y}  =  (pci)"'. 
y  =  {(2ma  +  kr)^ + 

There  is  a  constant  missing  in  (7.2)  due  to  the  omission  of  the  term  (cr^+  t^)  in  (7.1).  Also, 
(7.2)  do  not  provide  the  actual  wave  speeds  since  W  of  (7.1)  has  been  divided  by  a  constant. 
These  do  not  affect  the  analyses  which  follow  because  we  will  be  concerned  with  the  relative 
wave  speeds,  not  the  absolute  wave  speeds.  The  eigenvector  r  =  (r,,r2)  which  satisfles  the 
eigenrelation 

(G-»?l)r  =  0  (7.3) 

is  tangential  to  the  wave  curve  and  hence,  by  (3.3), 

r2  _  dr  _  2a  _  —  (2ma  +  kr)  TV 
r,  da  (2ma  +  ^T)Ty  2a 

The  upper  sign  is  for  the  Cj  wave  curves  and  the  lower  sign  for  the  C3  wave  curves.  The  right 
hand  sides  of  the  equation  are  homogeneous  in  a  and  t  which  implies  that  wave  curves  are 
similar.  As  stated  earlier,  the  wave  curves  for  c,  and  03  are  orthogonal  to  each  other.  Unless 
stated  otherwise,  all  discussions  in  the  sequel  are  for  the  wave  curves.  The  results  of  course 
apply  to  the  C3  wave  curves  with  minor  modifications  if  necessary. 

Equations  (7.4)  can  be  integrated  in  a  closed  form  [21].  It  gives  a  family  of  wave  curves  in 
the  (a,  t)  space.  There  are  straight  lines  passing  through  the  origin  which  are  themselves  wave 
curves.  We  call  the  straight  lines  wave  curves  the  “wave  line”.  They  are  determined  by  seeking 
the  solution  of 

T  _  dr 
a  da' 


The  wave  lines  are  then  given  by 


and 


-  =  m  ±  +  1  +  k.  (7.5) 

r 


We  see  that  depending  on  whether  k  is  less  than,  equal  or  larger  than  -  (m^  +  1)  there  will  be 
one,  two  or  three  wave  lines  (Fig.  1).  There  is  an  exception  for  k  =  —  1  at  which  we  have  only 
two  wave  lines. 

Along  the  wave  curve,  i.e.,  along  the  direction  of  the  eigenvector  r,  the  wave  speed  c  may  be 
increasing  or  decreasing.  If  f/  denotes  the  differentiation  of  rj  in  the  direction  of  r.  we  have 

7?  =  (Vrj)r. 


From  (7.2)  we  see  that 

T]  >0,  if  c  decreases, 
f]  <0,  if  c  increases. 

Differentiation  of  the  eigenrelation  (7.3)  gives  us 

(G„  -  i7d.,)r,  +  (C„  -  ?jd„)ry  =  0, 


=  G„r,r,. 


and  hence 


(7.6) 
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Fig.  1.  Classification  according  to  the  number  of  wave  lines  (WL)  and  inflectioa  lines  (IL). 


If  we  let 


a  =  a,,  r  =  a2, 


in  (7.1),  the  potential  function  W  is  homogeneous  in  a^,  of  degree  three  and  has  the 
property 


3^W(o„a2) 

dOiOjOk 


r,r/*  =  6W(r,,  r^). 


(7.7) 


Noticing  that 

‘‘  dotdot  ’ 

■  ^W(a„  02) 

''  3a,5o,5a* 


(7.8) 


we  obtain  from  (7.6),  (7.8)  and  (7.7)  that 

6W(r,,r,) 
T}= - ^ 

rkrk 


(7.9) 


Thus  the  sign  of  fj  is  determined  by  the  sign  of  W(r,,  Tj)  [11].  In  particular,  the  “inflection" 
point  on  the  wave  curve  at  which  f|  =  0  is  determined  by 

W(r„r2)=0.  (7.10) 


Since  the  wave  curves  are  similar,  the  inflection  points  on  the  wave  curves  form  a  straight  line 
passing  through  the  origin  which  is  called  the  “inflection  line”.  Equation  (7.10)  is  a  third  order 
polynomial  in  which  means  that  there  exists  at  least  one  inflection  line. 

We  now  investigate  the  number  of  inflection  lines.  Using  (7.1)  and  (7.4),,  we  write  (7.10)  as 

(l-ifc)A^  +  3A  +  2m=0.  (7.11) 


A 


dr 

da 


If  /c  =  1.  we  have  A  =  -  2m/3  and  A  =  “=,  the  latter  is  a  double  root.  If  k  #  1,  we  let 


A  = 


(7.12) 


We  then  have  one,  two  or  three  inflection  lines  if  A  is  positive,  zero  or  negative,  respectively. 
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To  have  three  inflection  lines  we  must  have 


or 


m~  + 


1 

\-k 


<0. 


From  the  above  discussion  we  classify  the  geometry  of  wave  curves  into  five  cases  according 
to  the  number  of  wave  lines  (WL)  and  inflection  lines  (IL)  (Fig.  1). 

Case  1:  k<  —(l  +  m~). 

Case  2;  -(1 +  m')sjk<0. 

Case  3  ;  0</c  £  1. 

Case  4:  l<A:<l+m~^ 

Case  5;  1  +  m~~  <  k. 

Case  2  and  Case  3  have  the  same  number  of  wave  lines  and  inflection  lines  but  they  have  quite 
different  wave  curve  geometry.  The  special  situation  in  which  m  =  0  has  been  studied  in  [8,9] 
where  Case  5  does  not  exist.  The  classification  given  by  Schaeffer  and  Shearer  [11]  is  based  on 
the  potential  function 


W(o,  T)  =  jcr^  +  fitrr  +  or*, 

in  which  a  and  b  are  the  parameters.  This  is  different  from  the  W  in  (7.1)  but  can  be  shown  to 
be  equivalent  if  we  apply  the  invariance  propeny  (ii).  In  the  {a,b)  plane,  Schaeffer  and 
Schearer  classify  the  wave  curves  into  four  regions.  It  can  be  shown  that  Case  1.  2,  3,  4  here 
corresponds,  respectively,  to  Region  4,  3,  2,  1,  in  [11).  Case  5  here  also  belongs  to  Region  2. 
The  boundaries  which  separate  the  five  cases  presented  here  appear  to  have  a  simpler  algebraic 
expression. 

Before  we  investigate  in  detail  each  case,  we  list  below  some  general  properties  of  the  wave 
curves  which  apply  to  at  least  two  or  more  of  the  five  cases.  At  the  point  where  the  wave  curve 
intersects  the  o-axis,  the  slope  of  the  wave  curve  is  obtained  from  (7.4)  as  (Fig.  2) 


—  =  -(Vl  +  m^  +  m),  at  a>0,  r  =  0. 
do 

=  (Vl  +  —  m),  at  cr<0,  r  =  0. 

This  is  independent  of  k  and  therefore  applies  to  all  five  cases.  We  remind  the  readers  that  all 


Fig.  2.  The  wave  curves,  wave  line  and  inflection  line  for  Case  1. 
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discussions  are  for  the  c,  wave  curves  which  are  represented  by  bold  curves  in  the  figures.  We 
see  that  the  absolute  value  of  the  slope  at  a  >  0  is  larger  than  that  at  a  <  0.  Moreover,  the 
product  of  the  two  slopes  is  -1.  which  means  that  the  two  slopes  are  orthogonal  to  each  other. 
Next,  the  wave  lines  given  by  (7.5)  apply  to  k  3:  -  (1  +  m^),  i.e.  to  all  cases  except  Case  1.  Let 
0,  and  62  be  the  angles  the  wave  lines  make  with  the  r-axis.  We  have 

tan  0,  =  Vm‘  +  1  +  fc  +  m  >  0,  for  all  k, 

tan  02  =  Vm“  +  1  +  k  —  m  ^  0,  for  k  $  -  1. 

It  is  readily  shown  that  +  d2<:z/2  or  >rr/2  if  A:  <0  or  >0.  Also,  0,  >  jr/4  except  when 

—  k>2m<2  and  62  >  rr/4  or  <n/4  when  k  >  2m  or  <2m.  Finally,  for  all  cases  except  Case  4, 

there  is  only  one  inflection  line.  The  slope  A  of  the  wave  curve  at  the  point  where  the  wave 
curve  intersects  the  inflection  line  is  the  only  real  root  of  (7.11)  and  is  given  by 

ASO,  if  kSl. 


(7.14) 


(7.13) 


in  which  A  is  defined  in  (7.12).  If  we  solve  (7.4)  for  t/o  in  terms  of  dr! da  (which  is  A),  we 


obtain 


T  l-A^-2mA  l  +  2A^4-(l-k)A'‘ 
a  kA  kA 


(7.15) 


in  which  the  second  equality  follows  from  (7.11).  With  A  provided  in  (7.14).  (7.15)  furnishes 
the  slope  of  the  only  inflection  line  for  Case  1,  2,  3,  and  5.  We  now  discuss  each  case 
separately. 


Case  1:  k<  -(1  +m^) 

In  Fig.  2,  the  bold  (or  thin)  lines  are  the  c,  (or  C3)  wave  curves  and  the  bold  (or  thin)  dotted 
line  is  the  inflection  line  for  Ci  (or  C3).  We  show  in  the  figure  only  typical  wave  curves.  Since 
wave  curves  are  similar,  other  wave  curves  can  be  obtained  by  enlarging  or  diminishing  the 
wave  curves  shown  in  the  figure.  The  r-axis  is  a  wave  line.  If  a  is  the  angle  the  inflection  line 
makes  with  the  positive  r-axis,  we  have  from  (7.15), 


tan  a  = 


_ ^ _ 

1  +  2A^  (1  -  k)A*  ’ 


(7.16) 


where  A  is  given  in  (7.14).  A  is  the  slope  of  the  wave  curve  at  the  intersection  of  the  inflection 
line.  It  is  also  the  slope  of  the  line  AB  on  which  W{o,  r)  =  0. 

The  arrows  on  the  wave  curves  and  the  wave  line  represent  the  direction  along  which  the 
wave  speed  c  decreases,  i.e.  17  >0.  This  is  determined  by  applying  (7.9).  For  instance,  to 
determine  the  direction  of  fj  >  0  at  point  P,  we  draw  from  the  origin  towards  the  region  W  >  0 
a  line  parallel  to  the  tangent  to  the  wave  curve  at  P.  This  determines  the  direction  of  the  arrow. 


Case  2:  -  (1  +  m^) ^k  sO 

There  is  one  inflection  line  and  (7.16)  which  gives  the  angle  a  applies  here  also.  Fig.  3.  There 
are  three  wave  lines,  one  of  which  is  the  r-axis  while  the  other  two  make  an  angle  0,  and  0- 
with  the  negative  r-axis.  Case  2  can  be  divided  into  Case  2A  and  Case  2B  depending  on 
whether  0^  <  0  or  >0  which  is  equivalent  to  Ac  <  - 1  or  >  -  1 .  Each  case  can  be  further 
subdivided  according  to  the  relative  magnitudes  of  a,  0,  and  0;.  We  have  (Fig.  4) 

Case  2 A,;  0<a<  -02  =  0i,  (Ac  =  -(l  +  m“)). 

Case  2A2:  0<a<  -02<0,. 

Case  2Av  0<  -03<q’<0,. 

Case  2B,;  0=  02<  a<  0,,  (Ac=-1). 

Case  2B2:  O<02<n'<0,. 

Case  2B,;  0  <  a  <  0,  <  0,. 

Case  2B4:  0  =  O'  <  02  <  01, 


(Ac  =  0). 


Cas*  2A, 
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Coa*  2A, 


‘  a 

I  / 


8,  -  -8,-  tan"  '  ( m) 


Com  2B, 


Cos*  2B2,  2B 


01  -  ton"'  (  2m  ) 

02  s  0 


Cos*  2B4 


Umbilic  lin« 


Fig.  3.  The  wave  curves,  wave  lines  and  inflection  line  for  Case  2. 

In  all  cases,  we  have  0,  <  n/2,  62  <  ^/4.  The  geometry  of  wave  curves,  wave  lines  and  the 
inflection  line  are  shown  in  Fig.  3.  Case  2Aj  is  not  shown  but  its  geometry  is  similar  to  that  of 
Case  2A2.  Case  2B4  is  a  degenerate  case  in  which  the  r-axis  is  an  umbilic  line.  It  is  also  a  wave 
line  for  c,  as  well  as  for  Cj. 

The  boundary  between  Case  2B2  and  Case  2B3  in  Fig.  4  is  characterised  by  the  relation 

a  =  02- 

With  a  given  by  (7.16)  and  (7.14)  and  02  by  (7.13)2,  this  boundary  can  be  plotted  numerically. 
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Fig.  4.  Further  classification  of  Case  2. 


It  can  be  shown  that  k—>—  1/2  as  m— *®.  An  alternate  way  in  finding  the  equation  for  a  =  62 
will  be  given  later.  [See  equation  (7.26).]  The  boundary  between  Case  2A2  and  2A3  is 
determined  by  the  condition 

ar  =  —62. 

Referring  to  Case  2A2  in  Fig.  3,  this  implies  that  the  wave  curve  and  the  inflection  line  are 
orthogonal  to  each  other  at  the  point  of  intersection  P.  With  dj  given  in  (7.13)2,  ibe  slope  A  of 

the  wave  curve  at  P  is  _ 

A  =  tan  02  =  Vm^  +  l+k  —  m. 

Substitution  of  A  into  (7.11)  leads  to  a  fourth  order  polynomial  in  k  and  m  which  can  be 
factorised  as 

[(2  +  k)^  +  4m^][(2  -  k)\k  +  1)  -  4m\k  -  1)]  =  0. 

The  first  factor  is  positive  and  non-zero.  We  therefore  have 

4m^  =  ^j^i2-kf.  (7.17) 

This  is  the  curve  denoted  by  or  =  -  ©2  F>g-  4. 


Case  3:  0<k£l 

Figure  5  shows  the  geometry  of  wave  curves,  wave  lines  and  the  only  inflection  line  which 
makes  an  angle  oto  with  the  negative  r-axis.  The  angle  Uq  is,  from  (7.15), 

—  kA 

~  l-t-2A^-l-(l-ik)A‘‘ ' 

and  A  is  given  in  (7.14).  It  can  be  shown  that 

0  <  OTo  <  02  <  <  ^/2. 

As  k— »0,  Uo-^O  and  the  r-axis  becomes  an  umbilic  line.  At  A:  =  1, 


tan  Oo  = 


6m 

9  -l-  8m^  ’ 


(7.18) 


and  c,  on  the  positive  r-axis  and  Cj  on  the  negative  r-axis  assume  the  same  constant  value.  The 
r-axis  is  therefore  a  shock  wave  curve. 
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Fig.  5.  The  wave  curves,  wave  lines  and  inflection  line  for  Case  3. 

Case  4:  1  <k  <  1 

This  is  the  most  complicated  of  all  five  cases.  It  has  three  wave  lines  as  well  as  three 
inflection  lines  (Fig.  6).  The  location  of  three  inflection  lines  are  obtained  by  the  slope  A  of  the 
wave  curve  at  the  point  of  intersection  with  the  inflection  line.  The  slope  A  is  a  root  of  (7.11). 
For  Case  4,  A  of  (7.12)  is  negative  and  (7.11)  has  three  real  roots.  Denoting  the  three  roots  by 
Ao,  A,,  A2,  we  have 

Ao  =  2(/c  -  1)“"^ cos(-  3^^)  , 

A,  =  2()t-l)-''^cos(^^^^). 

A2  =  2(A:-1)~‘^cos^j^, 


Fig.  6.  The  wave  curves,  wave  lines  and  inflection  lines  for  Case  4 
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If  tto,  Oi.  012  are  the  associated  angles  the  inflection  lines  make  with  the  r-axis.  Fig.  6,  we 
obtain  from  (7.15), 


-  kAj  ^  ^ 

tan  = - ; - T,  1  =  0,2, 

1 +  2Af +  (1 -k)A-‘ 


(7.19) 


tan  a,  = 


1  +2A?  +  (1  -*)At 


It  can  be  shown  that 


0<  ffo<  02<  01 
0  <  a,  <  (Xi , 

or,  +  02  <  Jt,  ai  +  0|  <  71. 


(7.20) 


Thus,  while  ao,  0|,  02  are  less  than  nil,  or,  and  a2  can  be  larger  than  nil. 

Depending  on  the  relative  magnitudes  of  a^,  or,,  0^2*  0i  and  02,  Case  4  can  be  subdivided  as 
shown  in  Fig.  7.  The  subdivisions  are  obtained  by  the  curves  in  the  {k,  m)  plane  which 
represent  or,  =  Oq,  oc\  =  62,  ■  ■  ■  ,  etc.  Each  subregion  is  identified  by  one  capital  letter  and  one 
lower  case  letter.  The  capital  letters  represent  the  following  conditions; 


A:  0  <  ar2  <  oto- 
B:  aQ<  a2<  6z- 
C;  02'''^2'^0i- 
D;  0i<a2<nl2. 

E:  nll<a2<n~dt. 


Fig.  7.  Further  classification  of  Case  4  and  Case  5. 
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The  lower  case  letters  represent  the  following  conditions: 

a:  0  <  a,  <  ao- 
b:  ao<  Ofj  <  62. 
c:  62<a;^<6i. 
d:  d,<  ai<A;'Z. 
e:  7r/2<  a,<jr~  6,. 
f:  ;r  —  0,  <  Or,  <  ;r  —  62- 

Therefore,  for  instance.  Case  4Be  is  governed  by  the  conditions 

aQ<a2<d2  and  :r/2<  at<  jz  —  6,. 

These  conditions  are  in  addition  to  (7.20)  which  govern  Case  4.  Case  4Cf  and  Case  4De  have 
two  more  subdivisions.  We  have 

a2<:z-  Oi  for  Case  4Cf,  and  4De,, 

a2>  Jt-  for  Case  4Cf2  and  4De2. 

The  equations  governing  the  boundaries  between  the  subclassifications  denoted  by  ai  =  61, 
02  =  ao,  ■  ■  •  ,  in  Fig.  7  have  the  following  expressions.  They  are 


k  +  1 

4m}  =  (2  kf-,  for  cti  -  0, 

k  -  1 

or  0:2  =  02. 

(7.21) 

^  ,  (V^-2)^  , 

4m  =  — p,  for  or,  =  jr/2  or 

02  =  nl2. 

(7.22) 

4m^  =  (4- k){k  -  1),  for  0:2  =  ao 

or  02  =  n  —  Ox, 

(7.23) 

,  j  (37  -  30/c)  T  (9  -  2k)\j9  -  Sk 
‘  m-i) 

6 

II 

.2 

(7.24) 

The  derivation  of  (7.21)  is  similar  to  that  for  a  =  -  62  in  Case  2.  It  is  therefore  not  surprising 
that  (7.21)i  is  identical  to  (7.17).  Equations  (7.22)-(7.24)  are  derived  as  follows.  We  substitute 
dr/daof  (7.4)  into  (7.11)  which  results  in  the  following  cubic  equation  for  rio,  the  slope  of  the 
inflection  line. 


+  -04  =  0, 


(7.25) 


where 


y  =  r/a, 


Uj  =  2mk^(k  -  1), 


a2  =  2k\k-l){Am^  +  l). 


03  =  6nik[{4  -  3^:)  -  4m\k  -  1)], 

a^  =  {k-  4y  +  8m\4  -U)-  I6m*(k  -  1) 

=  [(Vit  -  2)’  -  4m^(\/k  ~  OIKVa  +  2)^  +  4mHVk  +  1)]. 


If  or,  =  nl2  or  02  =  nl2,  y  =  0  is  a  root  of  (7.25).  This  means  that  04  =  0  which  leads  to  (7.22). 
If  0^2  =  OTo  ■  ^2  =  JT  -  a, ,  (7.25)  must  have  a  double  root  for  y.  The  condition  for  the  cubic 

equation  ^1  .lj)  to  have  a  double  root  can  be  shown  to  be  given  by  the  equation 

{k  - \)[mHk  -  1)  -  1  J(4m-  -{k-  1)(4  -  /fe)]*  =  0. 

The  vanishing  of  the  first  factor,  i.e.,  k  =  l.  is  a  degenerate  case  of  Case  4Aa  in  which 
O'!  =  a'2  =  0.  /c  =  I  is  also  the  boundary  between  Case  3  and  Case  4.  The  vanishing  of  the 
second  factor  gives  the  boundary  between  Case  4  and  Case  5  at  which  u,,  =  -  o , .  The 
vanishing  of  the  third  factor  leads  to  (7.23). 
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When  at  =  ofo.  >  =  ±  cot  »o  are  roots  of  (7.25).  This  means  that  (7.25)  can  be  faaoiized  as 

a^y^  +  aiy^-aiy-at 

=  (y^  -  cot^ao)(a,y  +  04  tan'ao), 


from  which  we  obtain 

2  ^2 

tan  ao  =  —  =  — . 

aj  04 


By  a  direct  substitution,  it  can  be  shown  that  the  second  equality  gives  (7.24). 

Finally,  if  we  denote  the  left  hand  side  of  (7.25)  by  /(y),  we  have 

/(cot02)  =  O.  for  a,  =  ©2.  (7.26) 

/(  — cot0i)  =  O,  for  a,  =  ;r  —  6,  or  a2  =  0i,  (7.27) 


in  which  6,  and  62  are  given  in  (7.13).  It  should  be  pointed  out  that  (7.26)  also  applies  to  the 
curve  denoted  by  a  =  63  in  Fig.  4  for  Case  2. 

Two  limiting  cases  of  Case  4  are  worth  mentioned.  In  the  limit  k  =  1,  we  have  a,  =  aj  =  0 
and  tto  is  given  in  (7.18).  This  is  a  special  case  of  Case  3.  At  the  other  limit  k  =  1  +  m~^,  we 
obtain 

rro=  =  7r  —  ^  tan"'(m“'), 


tan  a2  = 


2(1 +  m^) 


tan  0,  =  2m  +  m“'. 


This  belongs  to  Case  5  which  is  discussed  next. 

Case  5: 1 +m~^£k 

As  shown  in  Fig.  8,  there  is  only  one  inflection  line  which  makes  an  angle  aj  with  the 
positive  T-axis.  From  (7.15),  we  have 


tar.  a2  = 


_ ~  kA _ 

l  +  2A^  +  (l-/k)A^’ 


and  A  is  given  in  (7.14).  CAse  5  can  be  subdivided  into  5A,  5C,  5D  and  5E  as  shown  in  Fig.  7. 
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The  conditions  governing  the  subdivisions  are  as  follows. 

Case  5A:  O<a3<02. 

Case  5C:  02  <0^2 <  Pl¬ 
ease  5D:  0,  <  ar2  <  .t/2. 

Case  5E:  :tl2  <a-,<.T  -  d^. 

Thus  ar2  is  less  than  nil  except  for  Case  5E.  The  boundaries  02  =  02,  »2  =  0i .  and  02  =  nil 
between  the  subdivisions  are  provided  in  (7.21).  (7.27)  and  (7.22),  respectively. 


8.  CONCLUSIONS 

We  have  presented  here  a  refined  classification  of  2  x  2  non-strictly  hyperbolic  systems  near 
an  isolated  umbilic  point  for  a  general  potential  function.  The  classification  is  important  in 
solving  the  Riemann  problem  for  arbitrarily  prescribed  constant  initial  conditions  at  t  =  0,  x  >  0 
and  constant  boundary  conditions  at  jc  =  0,  /  >  0. 
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CHAPTER  15 


Growth  or  Decay  of  Shock  Waves  in  the  Generalized 
Goursat-Riemann  Problem* 

T.  C.  T.  Tingt 
Tankin  Wangt 


Abstract.  For  plane  waves  propagating  in  an  elastic  half-space  z  >  0  due  to 
an  impact  at  the  plane  z  =  0,  the  problem  is  that  of  a  generalized  Goursat-Riemann 
problem  in  which  the  initial  conditions  at  f  =  0  are  constant  but  the  boundary 
conditions  at  z  =  0  are  not.  Ifa(z,f)  is  the  traction  vector  at  any  plane  x  =  constant. 
*(0,t)  depends  on  t.  Assuming  that  «(0,f)  is  linear  in  t.  i.e.,  s(0,t)  =  -*-sf  t  where 

and  are  constant,  the  Goursat-Riemann  problem  corresponds  to  the  special 
case  in  which  vanishes.  We  study  the  generalized  Goursat-Riemann  problem 
with  non-zero  Aj.  We  investigate  how  the  solution  to  the  Goursat-Riemann  problem 
is  affected  by  the  non-zero  a“.  In  particular,  we  examine  whether  the  shock  waves 
generated  in  the  Goursat-Riemann  problem  grow  or  decay  due  to  the  non-zero 

1.  Introduction.  The  equations  govermng  the  plane  waves  in  an  isotropic 
elastic  solid,  with  certain  restrictions,  can  be  written  as  a  system  of  4  x4  conservation 
laws  [1,2], 


U, -f{U),  =0.  (1) 

U  =  ((7,T,U,v). 
f(U)  =  (pu,pv.t,-)). 


"The  work  presented  here  is  supported  by  the  V.  S.  Air  Force  Office  of  Scientific  Research  under  Contract  AFSOR- 
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t  Department  of  Civil  Engineering,  Mechanics  and  Metallurgy,  University  of  Illinois  at  Chicago. 
Box  4348.  Chicago.  IL  60680. 
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In  the  above,  cr.r  are  the  normal  and  shear  stresses,  £,7  are  the  normal  and  shear 
streuns,  u,v  are  the  normal  and  transverse  velocities  and  p  is  the  mass  density.  For 
hyperelastic  solids  there  exists  a  potential  function,  or  the  complementary  strain 
energy  W{(t,t)  such  that  I3] 


=  W^,  7  =  Wr, 


(2) 


in  which  the  subscript  denotes  differentiation.  Introducing  the  2-vectors 


U  =  (u,v),  p  =  (£,7). 


(3) 


equation  (1)  can  be  written  for  continuous  solutions  as 


U,  -  AUt  =0,  U  = 


(4) 


A  = 


where  I  is  a  unit  matrix  and 


G  = 


0  pi 
G  0 


"err  rrrr 


(5) 


The  charactenstic  wave  speed  c  is  obtained  from  14] 


(cA  + 1)^  =  0,  i  = 


L-r/(^c)J 


(6) 


which  can  be  reduced  to 


{G-(pc2)-‘l}r=0. 


(7) 


Assuming  that  G  is  positive  definite,  there  are  two  pairs  of  charactenstic  wave 
speeds  ^C] .  :rc3  and  we  let 


Cl  ^  cz  >  0. 
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Had  we  considered  fully  general  plane  waves,  ( 1 )  would  be  a  system  of  6  x  6  conser¬ 
vation  laws  and  there  would  be  an  additional  pair  of  characteristic  wave  speed  ±C2 
which  can  be  shown  to  be  hneairly  degenerate  [21. 

It  should  be  pointed  out  that  (7)  is.  mathematically,  the  characteristic  equation 
with  characteristic  wave  speed  for  the  following  2x2  conservation  laws 

+  f(*)i  =0,  'I 

(8) 

s  =  (^r,r),  f(a)  =  (W^„,ir.).  J 

Equation  (8)  is  a  familiar  equation  in  fluid  dynamics  and  hence  the  results  obtained 
here  can  be  applied  to  problems  in  fluid  dynamics.  Many  studies  have  been  made 
on  (8)  which  is  in  general  not  strictly  hyperbolic  [5-9]. 

In  this  paper  we  consider  the  following  initial  and  boundary  conditions  for  the 
generalized  Goursat-Riemann  problem 

j(r,0)=#\  I  >  0,  I 

(9) 

t>0,  J 

where  are  constant.  The  superscripts  ”b”  and  "a”  stand  for  "before”  and 

"after”,  respectively.  The  constants  we  2-vectors  defined  in  (3)i.  The 

problem  reduces  to  the  Goursat-Riemann  problem  when  =  0. 

When  *1  =  0,  U(i,t)  depends  on  x/t  only  and  is  best  expressed  in  the  form  of 
wave  curves  in  the  U  space.  The  wave  curve  consists  of  a  sequence  of  simple  wave 
curves  and/or  shock  wave  curves  |10  .  If  the  boundary  conditions  are  prescribed  in 
terms  of  stress  s,  it  suffices  to  consider  wave  curves  in  the  sub-space  s  1.2].  This  is 
not  possible  if  the  boundary  conditions  are  prescribed  in  terms  of  velocity  [11.121. 
The  simple  waves  are  ’  centered  simple  waves”  because,  in  the  x-t  plane,  they  are 
represented  by  a  family  of  straight  hues  from  the  origin  x=t=0.  Along  each  straight 
line  U  is  a  constant.  The  shock  wave  in  the  x-t  plane  is  represented  by  a  double 
line  which  is  a  straight  hne  because  the  shock  wave  speed  is  a  constant.  Between 
the  simple  waves  and/or  the  shock  waves  we  have  constant  state  regions.  When 
sf  ^  0.  we  no  longer  have  simple  waves  and  constant  state  regions.  Moreover,  the 
shock  wave  speed  is  not  a  constsmt  and  may  grow  or  decay  as  time  increases.  We 
will  investigate  the  perturbation  of  the  solution  to  the  generalized  Goursat-Riemann 
problem  from  the  Goursat-Riemann  problem  in  the  region  near  x  =  0  and  t  =  0.  In 
particular,  we  exaxmne  whether  the  shock  waves  grow  or  decay  due  to  the  non  ^ero 
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jj .  The  perturbation  is  with  respect  to  small  x  and  t,  not  with  so  that  the  value 
ofa“  need  not  be  small. 

2.  Perturbation  in  the  constant  state  and  simple  wave  regions. 
For  the  generalized  Goursat-Riemann  problem  with  initi2d  and  boundary  conditions 
prescribed  in  (9)  we  consider  the  solution  U(i,t)  in  the  form  [13-15] 


A  =  xjt. 


where  ,  are  functions  of  A  only.  Clearly,  vanish  identi¬ 

cally  for  the  Goursat-Riemann  problem.  Differentiation  of  (10)  with  x  and  t  results 


Since  A  in  (4)  depends  on  U,  we  also  expand  A  as 

A(U)  =  A“(A)  +  tA'(A)  ^  +  •  • 

A°(A)  =  A(U“)  1 


aHA)  =  •VA° 


The  operator  V  is  the  gradient  with  respect  to  the  components  of  U.  With  (11) 
and  (12),  (4)]  leads  to 


(AA°  +  I)— U°=0, 

d\ 


( AA°  +  I)  =  A°U’  -  AA^ 

d\  dX 


We  will  study  the  first  order  perturbation,  namely,  the  solution  of  U' 
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Let  c°  be  a  root  of  the  determinant 


cOA^c^j  +  Ili  =  0. 


The  solution  of  U®,U^  depends  on  whether  A  equals  to  any  of  the  four  roots  c° 
or  not.  We  call  the  region  "regular"  if  A  c°  tind  "singular”  if  A  =  c°  For  the 
Goursat-Riemann  problem  the  regions  A  ^  c°  and  A  =  c°  correspond,  respectively, 
to  the  constant  state  region  and  the  simple  wave  region. 

2.1  Regular  region.  When  A  ^  c°,  (14)  and  (15)  yield 


=  constant, 

u'  =  {(A“)-'  +  ^i}q. 


where  q  is  an  arbitrary  constant.  Equation  (11)2  ^hen  gives 

Utii=o=(AVq.  (17) 

and  (16)  and  (17)  provide  the  relation 

U' =(AA°  +  I)U.l.=o  (18) 

which  allows  us  to  convert  Utit=o  to  and  vice  versa.  We  see  that  while  U°.  Utlt=o 
are  constant.  U'  is  linear  in  A. 

2.2  Singular  region.  When 


A  =  c«, 


equations  (14)  and  (6)i  imply  that 


where  is  a  proportionality  factor.  In  particular. 


(20) 
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The  proportionality  factor  (,’  is  determined  by  differentiating  (19)  with  A  and  using 
(20).  We  have 


1  =  -Cr^-Vc® 


(21) 


and  (20)  is  rewritten  as 


=  {r°  •  r°.  (22) 

dX 

It  is  clear  that  V  here  is  the  gradient  with  respect  to  the  components  of  s  because  c 
depends  on  not  on  u.  Integration  of  (20)  or  (22)  furnishes  the  simple  wave  curves 
in  the  stress  space  1 1.2]. 

Inserting  (13)2  in  (15)  we  obtain  the  differential  equation  for  U'  for  the  singular 
region  as 


(AA°  ^  I)  =  A“U‘  -  a  •  (VA°)4-U°  (23) 

dX  dX 

which  is  bnear  in  Hence  if  is  a  solution,  so  is  kU^  where  k  is  an  arbitrary 
constant.  Despite  the  fact  that  the  coefficient  of  dU^  /dX  in  (23)  is  singular,  it  is 
shown  in  [13]  that  (23)  can  be  solved  for  dV^/dX  and  that  the  solution  for  U*  is 
unique. 

3.  Discontinuity  of  U‘  between  the  singular  and  regular  regions. 
We  consider  U(z.t)  on  a  curve  F  in  the  x-t  plane  which  is  given  by 

,  1  2 
If  —  Af  f  ~  ^  'tr  ^ 

where  Aj,  and  are  constant.  The  curvature  of  i-  at  t  =  0  is  k^.  Hence 

A  —  If.  /t  =  Ap  —  Kpt  —  . . . 

and.  from  (10).  U(x.<)  on  F  is 

UlZp.fi  =  U‘’(Ap)^(  {U'  -  -O(f-). 

If  F  is  a  boundary  between  a  regular  region  and  a  singular  region  and  if  F  is  not  a 
shock  wave.  U(z.f  I  and  U°(Ap  |  are  continuous  across  F.  We  then  have 
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} re|ui»r  region  —  {U  o^r  }»ingular  region 


This  shows  that  is  in  generzJ  discontinuous  across  the  boundary  between  a 
regular  region  and  a  singular  region  even  when  U(z,t)  is  continuous. 

4.  Perturbation  of  the  shock  wave  solution.  The  Rankine-Hugoniont 
equation  for  (1)  is 

[Uj^  Vlf(U)l  =0  (25) 

in  which  I”  is  the  shock  wave  speed  and  the  square  brackets  denote  the  jump  across 
the  shock  wave.  Employing  the  2- vectors  (3),  (25)  can  be  written  as 


[si  f  pV'lu:  =  0 
iu;  T'’[p;  =0. 


Let  Xv(t)  be  the  locus  of  the  shock  wave  in  the  x-t  plane  and  be  given  by 

Z„  =  V  <  +  -! - 

where  V',  a,,  are  .  respectively,  the  initial  velocity  and  the  initial  acceleration  of  the 
shock  wave  at  t  =  0  The  shock  wave  speed  V‘  at  anytime  t  is 

V'  =  V  ^  (27) 


Equation  (10 1  written  for  u  and  s  are 


u(a:,t)  =  u°  -i-  tu^  -  y  u^ 
*(z,t)  =s°  +  ts^  ~jS‘ 


With  the  notations  (3)3,  (2)  can  be  wntten  as 
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and  the  expression  of  p  in  the  form  (10)  leads  to 


p(«l  =  pO  _  ^pI  ^ 


p'^  =  p-  =  GV, 

„2  _  ^0  2  1  1 


In  the  above,  G°  =  G(«°).  Substituting  (27)  -(29)  into  (26),  we  have 


js®;  —  pV’Iu®'  =  0 


fu°:  -  rip°;  =0 


for  the  zero  order  and 


[«' ■  -  pV'iu’'  -  Pfv,.(u°i  =  0 
(U^:^  V(GV;  +  k„(p°]  =0 


for  the  first  order. 

From  (30)  we  obtain 


and  the  use  of  these  results  in  (31)  leads  to 


pViu-  =  -(*•;  -  (,ijy)[«0j, 

iu^;  =  -v[gV:-(n„/>v^^)[/]. 


Elimination  of  lu']  between  the  two  equations  in  (33)  yields 


/iI'^^GV  -  V-  =  -2(/t^/K)f*°J. 


Let  s~  ,s  be  .  respectively,  the  value  of  s  in  front  of  and  behind  the  shock 
wave.  When  s~  is  given.  (32;;  provides  a  solution  for  with  V  as  the  parameter. 
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As  V  varies.  s~  traces  out  a  curve  in  the  stress  space  which  is  the  shock  wave  curve. 
Fig.  1.  Equation  (32)i  then  the  jamp  \i  u°  while  (34)  and  (33)i  or  (33)2  give 
the  jump  in  and  in  terms  of  Kv 

5.  Growth  or  decay  of  the  shock  wave.  When  the  initial  acceleration  k^. 
of  the  shock  wave  is  non-zero,  the  shock  wave  grows  or  decays  cis  it  propagates.  If  the 
state  ahead  of  the  shock  wave  is  at  rest  or  at  a  uniform  state,  the  stress  amplitude 
and  velocity  amplitude  in  general  increase  or  decrease  according  to  whether  the 
acceleration  /«,.  is  positive  or  negative.  If  the  state  zdiead  of  the  shock  wave  is  not 
a  uniform  state,  the  positivity  of  may  not  be  accompanied  by  the  growth  of  the 
stress  amplitude  and/or  the  velociy  amplitude.  This  was  pointed  out  in  [16]  for  one¬ 
dimensional  shock  waves  and  in  [17-19]  for  three-dimensional  shock  waves.  Thus  the 
definition  of  growth  or  decay  of  a  shock  wave  depends  on  whether  one  is  looking  at 
the  stress  amplitude,  the  velocity  amplitude  or  the  acceleration  k^. 

In  this  paper  we  define  the  growth  or  decay  of  shock  waves  according  to  whether 
/Cf  is  positive  or  negative. 

6.  The  solution  to  the  generalized  Goursat-Riemann  problem.  The 
perturbation  presented  in  Sections  2,  4  and  the  discontinuity  condition  of  be¬ 
tween  the  singular  and  regular  regions  derived  in  Section  3  are  the  basis  for  the  con¬ 
structing  the  solution  to  the  generalized  Goursat-Riemann  problem  prescribed 
in  (9).  As  an  illustration,  consider  the  special  case  in  which  in  (9)  are  located 
on  a  shock  wave  curve  associated  with  C]  which  is  denoted  by  Vj  in  Fig.  1.  If  a“  =  0. 
the  solution  in  the  x-t  plane  consists  of  a  single  shock  wave  of  wave  speed  Vi ,  Fig. 
2.  Since  »(  —  0.  we  would  have  the  wave  pattern  shown  in  Fig,  3  in  which  there  is  a 
non-constant  shock  wave  Vj*  followed  by  an  acceleration  wave  C3.  The  acceleration 
wave  C]  does  not  appear  because,  by  Lax  stability  conditions  [20  .  the  C]  wave  speed 
behind  the  shock  wave  is  larger  than  the  shock  wave  speed  I'j’.  The  region  between 
the  ly  shock  wave  and  the  C3  acceleration  wave  denoted  by  M  is  a  regular  region 
which  is  not  a  constant  state.  We  will  determinate  the  imtial  acceleration  of  the 
I  j’  shock  and  sj"  which  is  the  value  of  in  the  region  M  as  t  —  0.  In  the  sequel,  the 
superscripts  a.b.m.---  denote  the  value  of  the  quantity  concerned  in  the  regions 
A.  B.  M.  ■  ■  ■ .  respectively,  as  t  — •  0. 

The  first  step  is  to  find  U'  from  U’  across  the  ohock  wave  I'j".  Applying 
(34)  to  the  shock  wave  I  j'  at  t  =  0  and  noticing  that  tne  reeion  B  is  a  constant 
state  which  means  that  U’  =  0,  we  have 

pV?G‘'a^~  =  -2(n,.  Tj)  V-' .  (35 1 
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The  regions  M  and  A  are  both  regular  regions  and  hence  which  leads  to 

G°  =  G”*.  Let 


=7ir]+7rr5,  (36) 

(a°l=a“-/  =  -(7?r?  +  7?rS),  (37) 

in  which  7f  of  i7ii73  ^^e  constant  and  Fj.rj  are  the  right  eigenvectors  of  (7)  as¬ 
sociated  with  c  =  cf,C3,  i.e., 


Mc?)2G‘‘-I}r“=0,  I 
Mc?)2G‘‘-I}r5=0.  j 


(38) 


In  the  above  it  is  understood  that  =  c®  and  =  Cj.  The  direction  of  r"  and 
Fj  are  taken  to  be  in  the  direction  of  decreasing  ci  and  C3,  respectively.  With 
given  in  (9),  7°  and  7°  are  determined  from  (37).  For  the  V)  shock  wave  and  the 
rj,r3  shown  in  Fig.  1,  we  have 


7?  >  0.  (39) 

Using  (36)  -(38),  (35)  can  be  written  as 

{(V,/ct)'  -  l}7rr?  ^  {(Vi/c5)'  -  l}73-r5 
=  2(^„/Vi)(7?f!  ^7?r3)- 


Hence  7j  and  73  are  given  by 

{(Vi/c»)2  -  l}7r  =2(«»/Vi)7?  ' 
{(Vi/c5)2-1}73-  =2(«./V3)73°  , 


(40) 


with  Kf  as  an  unknown  parameter.  We  have  thus  obtained  of  (36).  As  to 
using  (36)  and  (37).  we  write  (33)i  as 


^73'r?)  -  l«r/U,)(7?r?  *73°r?). 
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We  are  interested  in  /3V’]U^  —  s'"  ,  not  in  pV’iU^  .  Therefore  we  employ  (36)  and 

(40)  to  obtain 


pV'iu’ 


=  -(V^^—  WJ  -  (Vi/cS)i-l 


The  next  step  is  to  convert  to  U7‘  using  (18)  which  is  written  in  full  2is 


=s^^pViu'^, 


=  Vi  GV  +  u? 


a^m  ,  .,m 


Eliminating  u"  leads  to 


pViu^"-s^  =  (pVi^G"  -  I)^^-  (42) 

To  determine  *7^  expHcitly,  we  let 

*t  =1\  ^1+^3  ^3'  *«  -  7l  fl  +  73  f3> 

where  'r7‘,7r>7S‘i73  constant.  With  prescribed  in  (9),  7f  and  7J  are  de¬ 
termined  from  (43)2  The  boundary  between  the  region  A  and  region  M  is  a  C3 
charactenstic.  This  means  that  the  discontinuity  in«7  and  s'P  must  be  proportional 
to  r®.  Hence 


7i  -7i- 


(44) 


Equation  (42)  with  the  aids  of  (38)  and  (43)]  has  the  expression 


pViu'  -a"'  -  {(Vi/c7)*  -  l}7r'-?  +  )'  -  l}73"*r? 

and  from  (41 )  we  obtain 


a_  (V^l/ct)^  +  3 

-  {(Vi/c?)2-ip  Hi  VI  ’ 

m  _  (I'l/Cf)*  ~  3  \  0 

-  {(Vi/C“)‘-1}-  W]  V3- 


(451 
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Equation  (45)i  provides  the  initiad  acceleration  itt,  of  the  shock  wave  Fj*.  Since 
is  positive  by  (39),  and  7“  ha\'  the  opposite  signs.  In  other  words,  the  shock 
wave  Vj’  grows  when  7“  <  0  and  decays  when  7“  >  0.  For  the  shown  in  Fig.  1. 
7^  <  0  and  the  shock  wave  Vj*  grows  as  indicated  by  the  positive  curvature  of  Vj* 
in  Fig.  3.  The  Vcdue  of  73  has  no  effects  on  the  growth  or  decay  of  shock  wave  IjV 

If  -y®  =  0,Kr  =0  and  j”  =  0.  The  wave  pattern  in  the  x-t  plane  is  then 
given  in  Fig.  4  in  which  the  region  M  is  a  constant  state  and  the  region  A  has  a 
non- centered  simple  wave  solution. 

7.  Concluding  remarks.  The  example  given  in  the  last  section  illustrates 
how  one  can  determine  the  growth  or  decay  of  a  shock  wave  in  the  generalized 
Goursat-Riemann  problem.  The  example  was  one  of  the  simplest  cases  in  which  the 
two  stress  states  and  j**  are  on  a  shock  wave  curve  Vj .  Other  simple  c<ises  have  the 
two  stress  states  on  a  shock  wave  curve  V3  or  on  a  simple  wave  curve  cj  or  C3.  For 
the  latter,  the  discontinuity  condition  (24)  for  will  be  needed.  More  complicated 
cases  would  have  the  two  stress  states  at  arbitrary  locations  in  the  stress  space  in 
which  can  be  reached  from  through  a  series  of  simple  wave  curves  and/or  shock 
wave  curves.  More  studies  on  the  subject  is  being  undertaken  and  results  will  be 
reported  in  the  future. 
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